The heavy quark's self energy from moving NRQCD on the lattice 
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We present a calculation of the heavy quark's self energy in moving NRQCD to one-loop in perturbation 
theory. Results for the energy shift and external momentum renormalisation are discussed and compared with 
non-perturbative results. We show that the momentum renormalisation is small, which is the result of a remnant 
of re-parameterisation invariance on the lattice. 



1. INTRODUCTION 

The study of semi-leptonic B decays is cur- 
rently of great interest as it provides us with a 
test of the Standard Model. A key process is that 
of -B ^ TT semi-leptonic decay, however, because 
the B meson is much heavier than the final states, 
the recoil momenta of the tt meson can be large. 
This leads to both problems of small signal to 
noise ratios and significant discretisation errors 
if apir is large. Computations involving the sim- 
ulation of heavy quarks at low momentum can 
be dealt with using the NRQCD formalism PP, 
where the dynamics of interest are isolated by 
re-expressing the Lagrangian as an expansion in 
v/c. Although this framework for dealing with 
heavy quarks reduces the impact of discretisation 
errors, it is not designed to deal with large mo- 
menta. In this case, a suitable formalism to use 
is moving NRQCD (niNRQCD) pSll] which al- 
lows the momentum to be shared more evenly 
between the initial and final states. mNRQCD is 
developed in two stages, the first of which involves 
a choice of lattice frame such that the B meson 
is moving, thus reducing the magnitude of the 
pion's recoil momentum. Secondly, an effective 
field theory must be constructed for the heavy 
quark with non-zero velocity. In mNRQCD, the 
heavy quark's momentum is written as 



where ' 



- 7(1, v) is the 4- velocity of the b quark 
which is treated exactly, and is the small in- 



ternal momentum (of the heavy quark) inside the 
B meson which is discrctised on the lattice. 

As with the non-moving case, mNRQCD is 
constructed from a set of non-renormalisable in- 
teractions, each of which is accompanied by a 
coupling coefhcient which can be calculated per- 
turbatively. The heavy quark's self energy in 
NRQCD has been calculated to 0{as) Here 
we present a calculation of the one-loop renor- 
malisation parameters that appear in the mN- 
RQCD Lagrangian. In addition to the overall 
energy shift, mass and wavcfunction renormali- 
sation, we must also now calculate the renormali- 
sation of the external momentum. This renormal- 
isation is close to one because of the remnants of 
re-parameterisation invariance on the lattice. 

2. COMPUTING THE SELF ENERGY 

In this section, we present the Lagrangian used 
in this calculation and define the renormalisation 
parameters. 

2.1. Perturbation theory 

The Lagrangian is defined as 

ajC-NRQCB ^ ■iJjHx)ipix) - ip\x + at)llTpix) (1) 
where 

and 

i/„ = -^(A(^)-[v.V]^)-.v.V 
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In Eqn.^ n is a positive parameter that stabilises 
the high momentum modes in the heavy quark 
propagator. The mean-field parameter uq is used 
to improve the match between lattice and con- 
tinuum theories. The gluon fields are discretised 
using the Wilson plaquette gauge action. This 
Lagrangian includes all spin-independent correc- 
tions up to 0(1/M, 

The Feynman rules are calculated by expand- 
ing the U fields in terms of the coupling g. Follow- 
ing a Fourier transform to momentum space, the 
quark (gluon) propagator is defined as the inverse 
of the term that is bilinear in the quark (gluon) 
fields. Vertices of different order are selected ac- 
cording to their order in g. The rules are defined 
in the Feynman gauge and the infra-red diver- 
gence is regulated by introducing a small gluon 
mass. Note that in the limit — > 0, we recover 
the Feynman rules for the corresponding NRQCD 
action. 

2.2. Outline of the calculation 

The heavy quark's self energy S, represented 
diagrammatically in Fig. ^ along with tadpole 
counterterms from the uq factors, is defined by 
writing the inverse quark propagator G~^{k) in 
the form 

G'-i(fc) Q~Hk) - aE(fc) (2) 

where is the free propagator. The self energy 
can also be expanded for small A; |S] to give 



aT.{k) « Q!s ( fig ~ iko^li 

k2 
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(3) 



Combining the expansion of S with the quark 
propagator, to 0{as), gives 

G {k) ~ —ikf) — agilo + Oigikoili + 
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as^^^f^2 + v.k - a^v.krj^ H (4) 



Re-arranging this expression, we obtain 



k 



(5) 



Figure 1. Diagrams contributing to the self en- 
ergy. The broken lines represent temporal gluons 
and the curly lines represent spatial gluons. 

The renormalisation parameters are then defined 
as 



= 1 - as{flo + ill) 



(6) 



Zm = i + as{n2-ni) + as{n,~ni)v^-f^ (?) 

Zp = l-a,(a. -^^2) (8) 

where Mr = ZmM and Pr = ZpP (P = 7i\/v). 
The energy shift is defined as Eq = —agflo- We 
also note that since the tadpole counterterms are 
the same for ^2 and ily, these factors will cancel 
in Zp. In order to compute the renormalisation 
parameters, we need to isolate the coefficients, 
fin, of individual terms appearing in Eqn.|2| This 
is equivalent to taking derivatives of the diagrams 
with respect to an appropriate choice of momen- 
tum. 

The quark propagator is zero for t < 0, so by 
choosing the momentum flow appropriately, poles 
in the quark propagator can be avoided. The 
diagrams have internal loop momentum q and 
contain a pole in the gluon propagator in the go 
plane. The integrals are therefore evaluated an- 
alytically in the complex plane (z = e*''") which 
leads to a 3-D integral that can be evaluated nu- 
merically using the VEGAS package. The calcu- 
lations were performed at aM = 2.0 with several 
different values of the stabilisation parameter and 
velocity. These values match those used in the 
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Figure 2. The simulated energy (E'sim) and the 
physical binding energy ((E'sim — Eo)jr) for a 
heavy- light meson with aM = 2. 

numerical calculations [S] which were carried out 
on Wilson glue configurations at (3 ~ 5.7. We 
take as{l/a) = 0.35. The results are checked for 
invariance under Vj —Vj. 

3. PRELIMINARY RESULTS 

Preliminary results for the heavy-light binding 
energy Mkin— ^m-M = (i?sim— i?o)7r as a function 
of the bare velocity are presented in Fig. El The 
simulated energy [S] clearly exhibits velocity de- 
pendence as v/c approaches 1. After subtraction 
of the energy shift, the physical binding energy is 
independent of v, as it should be. 

Preliminary values for the renormalisation of 
the external momentum 7MV as a function of the 
bare velocity are presented in Fig.O Up to lattice 
artefacts, the Lagrangian has an invariance under 
arbitrary shifts of momentum between jMv and 
k, and this protects Zp from becoming very dif- 
ferent from 10. Zp from perturbation theory 
and simulations disagree at the level of system- 
atic errors in this calculation. 

4. CONCLUSIONS 

We have presented preliminary results for the 
energy shift and external momentum renormal- 



Figure 3. The renormalisation of the external mo- 
mentum as a function of bare velocity (v/c). Sim- 
ulation errors are statistical only. Perturbative 
errors are taken as 1 x a^w^. 

isation for a heavy quark in moving NRQCD. 
Comparison with numerical simulation is encour- 
aging. 

The full set of renormalisation parameters will 
be presented in a forthcoming publication |7] . Fu- 
ture work will incorporate higher order terms in 
the Lagrangian and include the renormalisation 
of current operators for B tt decay. 
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